EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 





CERN-PH-EP-2010-011 
2010/05/31 



CMS-MUO-10-001 



Measurement of the charge ratio of atmospheric muons 

with the CMS detector 



o 

(N 



The CMS Collaboration 



(N 



X 

43 



> 

(N 

m 
m 
in 

in 

o 
o 

> 

X 



Abstract 

We present a measurement of the ratio of positive to negative muon fluxes from cos- 
mic ray interactions in the atmosphere, using data collected by the CMS detector 
both at ground level and in the underground experimental cavern at the CERN LHC. 
Muons were detected in the momentum range from 5 GeV/c to 1 TeV/c. The surface 
flux ratio is measured to be 1.2766 ± 0.0032 (stat.) ± 0.0032 (syst.), independent of the 
muon momentum, below 100 GeV/c. This is the most precise measurement to date. At 
higher momenta the data are consistent with an increase of the charge ratio, in agree- 
ment with cosmic ray shower models and compatible with previous measurements 
by deep-underground experiments. 
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1 Introduction 

The muon charge ratio R is defined as the ratio of the number of positive- to negative-charge 
atmospheric muons arriving at the Earth's surface. These muons arise from showers produced 
in interactions of high-energy cosmic ray particles with air nuclei in the upper layers of the 
atmosphere. The magnitude and the momentum dependence of R are determined by the pro- 
duction and interaction cross sections of mesons (mainly pions and kaons), and by their decay 
lengths. As most cosmic rays and the nuclei with which they interact are positively charged, 
positive meson production is favoured, hence more positive muons are expected. Previous 
measurements from various experiments JTHSI showed the muon charge ratio to be constant 
up to a momentum of about 200 GeV/c, and then to increase at higher momenta, in agreement 
with the predicted rise in the fraction of muons from kaon decays. Measurements of the charge 
ratio can be used to constrain hadronic interaction models and to predict better the atmospheric 
neutrino flux. 

The Compact Muon Solenoid (CMS) |9| is one of the detectors installed at the Large Hadron 
Collider (LHC) IfTUl at CERN. The main goal of the CMS experiment is to search for signals of 
new physics in proton-proton collisions at centre-of-mass energies from 7 to 14 TeV |(TT| . 

Cosmic rays were used extensively to commission the CMS detector [T21IT3I . These data can 
also be used to perform measurements of physical quantities related to cosmic ray muons. This 
letter presents a measurement of the muon charge ratio using CMS data collected in two cosmic 
ray runs in the years 2006 and 2008. More details of the analyses can be found in |T4l[l5| . 



2 Experimental setup, data samples, and event simulation 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diame- 
ter, providing a field of 3.8 T Within the field volume are the silicon pixel and strip tracker ||T6|, 
the crystal electromagnetic calorimeter and the brass-scintillator hadron calorimeter. Muons 
are measured in gas-ionization detectors embedded in the steel return yokes [17j. In the barrel 
there is a Drift Tube (DT) system interspersed with Resistive Plate Chambers (RPCs), and in the 
endcaps there is a Cathode Strip Chamber (CSC) system, also interspersed with RPCs. In ad- 
dition to the barrel and endcap detectors, CMS has extensive forward calorimetry A detailed 
description of CMS can be found in H). 

The CMS detector is installed in an underground cavern, with the center of the detector 89 m 
below Earth's surface, and 420 m above sea level. The location is 46° 18.57' north latitude and 
6° 4.62' east longitude. The upper 50 m of the material above CMS consists of moraines, fol- 
lowed by 20 m of molasse rock. A large access shaft with a diameter of 20.5 m rises vertically to 
the surface, and is offset from the center of CMS by 14 m along the beam direction. It is covered 
by a movable concrete plate of 2.25 m thickness. Thus, depending on the point of impact on 
CMS, the total material traversed by close-to-vertical muons changes from approximately 6 to 
175 meters of water equivalent. 

The CMS experiment uses a right-handed coordinate system, with the origin at the nominal 
proton-proton collision point, the X axis pointing towards the center of the LHC ring, the y 
axis pointing upwards (perpendicular to the LHC plane), and the z axis pointing along the 
anticlockwise beam-direction, at geographic azimuth 280.8° (approximately west). The angle 
between the CMS y axis and the zenith direction is 0.8°. This small difference is neglected in 
the analysis, and the angle of the muons relative to the y axis is used to represent the zenith 
angle Z . 



3 Cosmic-muon reconstruction 



At the center of the detector, the magnetic field is parallel to the central axis of the solenoid, 
which is aligned with the z axis. Muon momenta are reconstructed by measuring the curvature 
of the muon trajectory projected on the xy plane, which yields the component of muon mo- 
mentum transverse to the z axis, pj = p sin 9, where 9 is the polar angle with respect to the z 
axis. This configuration is favourable for the reconstruction of atmospheric muons, providing a 
strong magnetic bending for muons traversing the detector, at any incident azimuthal angle (p 
around the z axis. Full tracking of muons is available in the polar angle range 10° < 9 < 170°. 

CMS collected cosmic ray data in several runs during the final years of detector construction 
and commissioning. Data from the Magnet Test and Cosmic Challenge in 2006 (MTCC) ffT2| 
and the Cosmic Run At Four Tesla in 2008 (CRAFT08) [13] are used in the analysis reported 
here. 

In August 2006 the CMS detector was pre-assembled on the surface before being lowered into 
the cavern. In this configuration no material above the detector was present, apart from the thin 
metal roof of the assembly hall. A small fraction of each of the subdetectors was instrumented 
and operating at the time. The details of the MTCC setup are described in |(T2l |14| . About 
25 million cosmic-muon events were recorded during the first phase of the MTCC with the 
magnet at a number of field values ranging from 3.67 to 4.00 T 

The CRAFT08 campaign was a sustained data-taking exercise in October and November 2008 
with the CMS detector fully assembled in its final underground position. The full detector, 
ready for collecting data from LHC, participated in the run, with the magnet at the nominal 
field of 3.8 T Approximately 270 million cosmic-muon events were recorded. 

Single cosmic muons are simulated using the Monte Carlo event generator CMSCGEN |[T8"IIT9"I , 
which makes use of parameterizations of the distributions of the muon energy and incidence 
angle based on the air shower program CORSIKA [20 1. The CMS detector response is simulated 
using the GEANT4 program [21j, which takes into account the effects of energy loss, multiple 
scattering, and showering in the detector. A map [19 j describing the various materials between 
the Earth's surface and the CMS detector is used to obtain the average expected energy loss 
of simulated muons as a function of their energy, impact point, and incidence direction at the 
surface. 



3 Cosmic-muon reconstruction 

Muon tracking in CMS can be performed with the all-silicon tracker at the heart of the detec- 
tor, and with either three or four stations of muon chambers installed outside the solenoid, 
sandwiched between steel layers serving both as hadron absorbers and as a return yoke for the 
magnetic field. 

Three types of muon-track reconstruction were designed for cosmic muons not originating 
from an LHC proton-proton collision [22J: a standalone-muon track includes only hits from 
the muon detectors; a tracker track includes only hits from the silicon tracker; and a global- 
muon track combines hits from the muon system and the silicon tracker in a combined track 
fit. For a cosmic muon that crosses the whole CMS detector, illustrated in Fig. [I] (top), each 
of the above types of tracks can be fitted separately in the top and bottom halves of CMS. 
Alternatively, a single track fit can be made including hits from the top and bottom halves of 
CMS. The direction of the muon is assumed to be downwards, and the muon charge is defined 
accordingly. 

The analysis based on 2006 MTCC data uses standalone muons. The reduced detector setup 
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Figure 1: Cosmic-ray muons crossing the CMS detector. The upper two pictures display 
muons from 2008 underground data, leaving signals in the muon system, tracking detectors 
and calorimeters. A standalone track (top left) and a pair of global half-tracks (top right) are 
shown. The bottom plot depicts a muon from 2006 surface data crossing the muon chambers 
at the bottom of CMS. 

used in the MTCC was just a fraction of the bottom half of the complete detector, depicted in 
Fig. FT] (bottom). Since the muons were measured only in one half of the detector, the momen- 
tum resolution is poorer than in the standalone-muon analysis using the complete detector. 
Having the detector on the surface, however, permitted the collection of a large number of 
low-momentum muons, down to a momentum of 5 GeV/c, allowing for a precise measurement 
of the charge ratio in the low-momentum range. 

Two analyses based on the 2008 CRAFT08 underground data are performed, one using stan- 
dalone muons and the other using global muons. The underground global-muon analysis 
(GLB) profits from the excellent momentum resolution and charge determination of global- 
muon tracks, but requires that the muon passes through the silicon tracker. The underground 
standalone-muon analysis (STA) profits from the larger acceptance of the muon chambers and 
yields approximately eight times as many muons as the global-muon analysis. In a standalone 
cosmic-muon fit spanning the whole diameter of the muon detector (Fig. IT), the momentum 
resolution is significantly improved compared to a standalone fit using only half the detector. 
The improvement varies from a factor of four at low momentum to more than a factor of ten 
for momenta above 100 GeV/c l|22|. 

The "maximum detectable momentum", p m dm/ defined as the momentum for which the cur- 
vature of a muon track is measured to be one standard deviation away from zero, is around 
200 GeV/c for standalone-muon tracks in one half of the detector, around 10 TeV/c for standalone- 



3 Cosmic-muon reconstruction 



muon tracks traversing the entire detector, and in excess of 20 TeV/c for global-muon tracks. 
The distribution of the transverse momentum (Pj CA ), calculated at the point of closest approach 
(PCA) to the nominal proton-proton collision point taking into account the energy loss in the 
detector, is depicted in Fig. [2] (a) for the muons selected in the global and standalone-muon 
underground analyses. 

The redundancy of the different tracking systems in the complete CMS detector allows the 
determination of the momentum resolution and rate of charge misassignment (the fraction of 
muons reconstructed with incorrect charge) directly in data. In the global-muon analysis, the 
half-difference of the track curvatures measured in the top half and the bottom half of the 
detector d.Q is used to measure the resolution of the half-sum Cj: 
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where px is the transverse momentum and q the charge sign of the muon. Both the core and the 
tails of the resolution distribution are well reproduced by the dc T estimator, as demonstrated 
for simulated events in Fig. [21(b). 
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Figure 2: (a) Normalized muon px distributions, for the global (closed circles) and standalone- 
muon analyses (open circles), at the PCA. Differences in the distributions are expected, as the 
global and standalone-track fits have different momentum resolutions and acceptances, (b) 
Comparison of the (q/pj) resolution estimate dc T (closed circles) with the true Cx resolution 
(hatched histogram), obtained from simulated global muons. 

In the underground standalone-muon analysis, an independently reconstructed tracker track 
is available in 40% of the selected events. The comparison of the charge and momentum mea- 
sured for the tracker track and for the standalone-muon track gives a measure of the tracking 
resolution both in data and in simulated events. 



All three analyses measure the charge ratio in events with a single cosmic ray muon, rejecting 
events with more than one muon detected. 

4 Event selection and analysis 
4.1 Analysis of surface data 

The cosmic-muon charge ratio was measured by CMS for the first time using MTCC data fT4 | . 
For this analysis, only the bottom sector in two (out of five) wheels of the barrel muon system 
(DT) is used. Selection accepts only muons triggered and reconstructed in a perfectly left-right 
symmetric fiducial volume with respect to the vertical axis, emphasized in Fig. IT] (bottom), 
ensuring a charge-symmetric acceptance. 

Around 15 million events were recorded in runs with DT triggers and a stable magnetic field 
above 3.67 T. About 330 000 events pass the fiducial-volume and track-quality selections. The 
measured muon charge ratio and its statistical uncertainty are displayed in Fig. [3] (a), as a func- 
tion of the measured muon momentum, before any correction due to detector effects is applied. 
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Figure 3: Uncorrected charge ratio, together with the statistical uncertainty, (a) From 2006 
MTCC data, as a function of the measured muon momentum, (b) For the global (closed circles) 
and standalone-muon analyses (open circles), as a function of the measured px at the PCA. 

The probability of charge misassignment is small for low-momentum muons. At high mo- 
menta, resolution effects increase the chance of charge misassignment thus lowering the mea- 
sured value of the charge ratio. Only muons with a measured momentum below p m dm = 
200 GeV/c are included in the analysis. 

4.2 Underground global-muon analysis 

The 2008 data were recorded using a single-muon trigger requiring the coincidence of muon 
hits in at least two muon detector layers. Triggers from the DT or RPC systems in the top or 
the bottom halves of the detector were accepted. The trigger efficiency is high for muons with 



4 Event selection and analysis 



sufficient momentum (a few GeV/c) to penetrate several layers of the steel return yoke 
The subsequent event selection is designed to ensure good track quality and high efficiency. 

The muon trajectory in each half of the detector is required to contain at least 20 (out of 44 
possible) hits in the DT system. Of these 20 hits, at least 3 hits are required to measure the 
longitudinal coordinate (z), ensuring a good measurement of the polar angle. The muon tra- 
jectory is required to contain no hits in the muon or tracker endcaps. The two halves, top and 
bottom, of each cosmic-muon trajectory are required to be reconstructed as two separate track 
segments in the silicon tracker, each containing at least 5 hits (out of 12 possible) in the tracker 
outer barrel system. A loose cut is applied to the normalized x 2 of each of the two global-muon 
fits and the polar angles are required to match within |A cot 9\ < 0.2, in order to suppress the 
small background from multi-muon cosmic shower events. The average transverse momen- 
tum of each muon, measured at the PCA, is required to be greater than 10 GeV/c in order to 
ensure that the muon is able to traverse the entire CMS detector. All selection requirements are 
applied to the top and bottom muon trajectories. 

While the main shaft of the CMS underground area is symmetric with respect to the yz plane, 
the two auxiliary access shafts are located at asymmetric positions with respect to this plane 
(cf. Section [2J. This causes the geometrical acceptance of the detector to be asymmetric for 
muons of different charges, since the CMS magnetic field is aligned with the z axis. To remove 
this effect, muon tracks that cross these auxiliary shafts are not considered in the analysis, nor 
are muons that cross the mirror images of those regions with respect to the x = plane. We 
refer to this requirement as "symmetric selection". 

About 245 000 muons are selected. The muon pj distribution is reported in Fig. [2] (a) for the se- 
lected muons. Figure 0(b) depicts the measured uncorrected charge ratio as a function of pj CA . 

4.3 Underground standalone-muon analysis 

In this analysis the particle trajectory is reconstructed using only the hits in the barrel muon 
system (DT and RPC). To select muon tracks that are fully contained in the barrel region, events 
with hits in the endcap CSCs are rejected. A single track is reconstructed using the information 
from both halves of the detector. Only one standalone muon per event is allowed. 

Muon tracks are required to have a transverse momentum, measured at the PCA, larger than 
10 GeV/c. At least 45 muon hits (out of 88 possible) are required to be associated with the track. 
The muon trajectory in the event is also reconstructed as two standalone-muon tracks, one in 
the upper and one in the lower half of the detector, with more than 20 hits (out of 44 possible) 
each. 

In order to ensure a good track-quality further selection criteria are applied to the tracks: the 
normalized x 2 of each reconstructed muon track must be less than 5, the impact parameter in 
the xy plane must be less than 100 cm, the track direction at PCA must be vertical within 42° 
in 9 and 60° in (p, and the track PCA must lie within the range \z\ < 600 cm. A "symmetric 
selection" is also applied as in the global-muon analysis. The number of muons selected is 
1.6 million. 

The analysis relies on the simulation to correct for charge misassignment and momentum reso- 
lution effects, using the data with both a standalone and a tracker track in the event to perform 
further corrections and estimate systematic uncertainties. From the comparison of tracks recon- 
structed both in the tracker and in the muon system, the probability of charge misassignment 
is known to be well below 1% for pj CA < 0.5 TeV/c, increasing up to about 1.5% in the highest 
momentum bin. The difference observed between data and simulation in the subsample of 



events that include a tracker track is taken into account to correct the charge misassignment 
and to assign the related systematic uncertainty, as explained in Section |6J 

The muon momentum scale and resolution are determined by comparing the transverse mo- 
mentum of the standalone-muon track to that of the associated tracker track, and are accurately 
modeled by the simulation. Therefore the momentum unfolding, which provides an estimate 
of the true momentum of the muon tracks from the measured momentum, can be based on 
the simulation. An uncertainty on the momentum resolution for all events, including those 
without a tracker track, is taken into account as a systematic uncertainty. The momentum scale 
in the tracker volume is set by the magnetic field, which is known to a precision better than 
0.1% |24], as confirmed by additional checks performed with early LHC data [25]. The uncor- 
rected muon charge ratio is shown in Fig. b^ (b) as a function of p j C 



„PCA 



5 Corrections for energy loss and resolution 

In order to express the charge ratio measurement as a function of the true momentum at the 
surface of the Earth, the measured momentum inside the CMS detector has to be corrected 
for energy lost between the surface of the Earth and the point of measurement. Furthermore, 
corrections need to be applied for migration of entries from bin to bin due to momentum reso- 
lution and for possible misassignment of the muon charge. 

5.1 Energy-loss correction 

In the MTCC analysis the measured muons are propagated back to the top of CMS, correcting 
for expected momentum loss and bending in the magnetic field. In addition, the effect of charge 
misassignment is estimated using simulated events, and a bin-by-bin correction is applied to 
the measured charge ratio. 

For the muons selected in the global and standalone-muon analyses of the 2008 underground 
data, the average expected energy loss depends strongly on the path followed through the 
Earth. The underground measurements are corrected for this effect by propagating the trajec- 
tory of individual muons back to the Earth's surface, using the same material model as in the 
simulation (cf. Sectional. Energy loss in matter is about 0.15% higher for }i + than for ji~ due 
to slightly larger ionization losses [7]. This difference is taken into account in the energy-loss 
correction, but affects the measured charge ratio by less than 0.3% over the entire momentum 
range. 

5.2 Unfolding the momentum spectrum 

In the underground data analyses, momentum resolution effects in the detector are corrected 
using an unfolding technique, applied to the charge-signed inverse momentum C = C]/p. 
In this procedure p represents the measured momentum extrapolated to the Earth's surface, 
where the correlation with the true muon momentum is highest. 

The momentum measured at the PCA is propagated first to the top of CMS, accounting for 
the magnetic field and the amount of material traversed, and then from the top of CMS to the 
surface of the Earth, following a straight line. The angular resolution of the detector is better 
than 5 mrad. Only muons with an estimated momentum above 30 GeV/c after this correction 
are kept in the analyses. 

Given a vector of true muon counts N true matrix inversion is used to compute the best estimator 



8 6 Systematic uncertainties 

Nj rue from the vector of observed muon counts ]Vf neasured : 

^measured = ^MijNf * , 

i 
Nf ue = ^M|T 1 Nj nea8ured . (2) 

The migration matrix element Mu is the probability that a muon with true C (C tme ) in bin j is 
observed with a measured C (c measured ) in bin /. M,y is an approximation of the exact migration 
matrix, and is constructed differently for the global and standalone-muon analyses. 

In the standalone-muon analysis the migration matrix estimator is extracted by comparing the 
true momentum to the reconstructed momentum in simulated events. 

In the global-muon analysis the approximate migration matrix is derived directly from the 
data. For each muon, the C values measured in the top and the bottom half of the detector 
are propagated individually to the Earth's surface. The estimated true C is then defined as 
(-true _ (Q op _|_ Cbottom)/2, and the measured values Qop and Cbottom are used to represent 
^measured They both have the desired property C measured = C true ±d c , where d c is the C reso- 
lution estimator, defined as dc = (C top — Q, ottom )/2. The matrix Mu is then populated using 
these estimated values, for all muons in the selected event sample. As the resolution estimator 
dc gives a good representation of the actual resolution of C tme (Fig. 0(b)), this procedure yields 
a good approximation of the true migration matrix Mjj. 

In both analyses, variations of the energy loss around the expected value are taken into ac- 
count in the unfolding procedure by applying an additional 10% Gaussian smearing of the 
energy-loss correction to the measured momentum when forming the migration matrix. This 
approximation is based on simulation studies using GEANT4. 

The muon counts N,- correspond to the bins of the histograms in which the corrected charge 
ratio results are presented. The bin boundaries were chosen such that the migration between 
bins is small. The values of the off-diagonal elements of the migration matrix are below 0.1 in 
the global-muon analysis and less than 0.2 in the standalone-muon analysis. 

The measurement of the charge ratio using 2006 data, corrected for energy loss in the detector 
and for charge misassignment, is depicted in Fig. B] as a function of the muon momentum, 
together with the statistical and systematic uncertainties. 

The measurements of the muon charge ratio in the global and standalone-muon analyses of 
2008 data are displayed in Fig. 15} as a function of the muon momentum. The "raw" result is 
based on the final alignment including the scale correction discussed in Section |6] The "cor- 
rected" results are based on the unfolding and, for the standalone-muon analysis, include an 
additional charge-misassignment correction. 

6 Systematic uncertainties 

Systematic uncertainties arise from reconstruction and instrumental effects that can affect dif- 
ferently the detection efficiency and momentum measurement of }i + and ]i~ . They are evalu- 
ated as a function of the muon momentum at the Earth's surface. 

The CMS magnetic field is known with high precision in the region inside the superconducting 
solenoid, and with less precision in the steel return yoke [24] . Systematic effects on the charge 
ratio due to the uncertainty on the magnetic field are less than 1%. 
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Figure 4: Charge ratio for the surface analysis, as a function of the muon momentum, corrected 
for energy loss in the detector and for charge misassignment, after propagating the muon track 
to the entry point in CMS. The thick error bars denote the statistical uncertainty and the thin 
error bars statistical and systematic uncertainties added in quadrature. 

A possible bias in the positive and negative muon rates detected underground, due to asymme- 
tries in detector acceptance and uncertainties in the material densities used in the material map 
(known within 5%), yields a non-negligible uncertainty on the charge ratio only in the lowest 
momentum bin. The additional effect of the selection cuts is generally small, well below 1%. 

The requirement of a muon trigger in the detector leads to a small difference in efficiency for 
positive and negative muons, below 1%, which is correlated between the two underground 
analyses. Both analyses estimate a possible systematic bias induced by the trigger by employ- 
ing a so-called tag-and-probe technique, using information from both halves of the detector 
and, in the case of the standalone-muon analysis, information from the independent DT or 
RPC muon triggers. 

In the global-muon analysis the effect of charge misassignment is small, ranging from less than 
0.01% at 10 GeV/c to about 1% at 500 GeV/c, and it is corrected by the unfolding procedure, 
using the data-driven resolution estimator defined in Eq. dlj. 

In the standalone-muon analysis the charge misassignment correction to the charge ratio, in- 
cluded in the unfolding matrix, is based on simulated events and tested in real data using the 
subsample of standalone muons with an associated tracker track. A higher rate of charge mis- 
assignment is observed in data than in simulation, with a maximum absolute discrepancy of 3% 
in the highest momentum bin. Since this discrepancy could not be attributed unambiguously 
to the standalone-muon tracks, a correction is applied equal to 50% of the full effect observed 
in data, with a systematic uncertainty equal to the correction itself. 

The precise alignment of all the tracking-detector components is crucial for accurate recon- 
struction of high-px muons, whose trajectories have only a small curvature in the detector. 
Cosmic muon tracks from the same 2008 data set used for this analysis are employed to per- 
form such an alignment of the silicon tracker and muon system [26. 27J. Possible effects from 
potential residual misalignment that could lead to momentum migrations and incorrect charge 
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Figure 5: Muon charge ratio as a function of the muon momentum at the Earth's surface for 
(a) global and (b) standalone-muon analyses. Open squares indicate the uncorrected ratio, 
including full alignment. Closed circles show the unfolded charge ratio with statistical errors 
only. The lines denote the statistical and systematic uncertainties added in quadrature. 

assignments are evaluated by studying various realistic missalignment scenarios in data and 
simulation. Only the two highest momentum bins are potentially affected by misalignment, as 
expected, yielding a bias in the charge ratio around 1% in the two highest-momentum bins for 
the global-muon analysis. For the standalone-muon analysis, the effect in the charge ratio is 
less than 1% up to 400 GeV/c, and around 4% in the highest-momentum bin. 

A global deformation of the detector could be missed during the alignment procedures (a so- 
called "x 2 -invariant" or "weak" mode 11281 ), and potentially affect the charge ratio. The most 
problematic deformation would be a mode which caused a constant offset in q / Pj CA , differ- 
ent from zero, affecting the momentum scale for cosmic muons of opposite charge in opposite 
directions. A two-parameter fit of the simulated q/p^ CA distribution to the data is performed 
using muons in the range py CA > 200 GeV/c, leaving the unknown charge ratio and the q/ Pj CA 
offset in the simulation to vary freely in the fit. An offset of 0.043 ± 0.022 c/TeV is found. The 
measured muon momenta are corrected for this offset and its uncertainty is included as an 
additional systematic uncertainty on R, fully correlated between the two underground mea- 
surements, of the order of 1% and 4% respectively in the two highest momentum bins. 

In the 2006 MTCC analysis, systematic uncertainties arise mainly from the finite precision of the 
detector alignment parameters, from the correction of the charge misassignment probability 
and from the slightly larger uncertainty (~5%) in the scale of the magnetic field in the steel 
return yoke. 

The total systematic uncertainties in the three analyses are summarized in TablellJ as a function 
of p at the Earth's surface. The systematic uncertainties have also been evaluated as a function 
of the vertical momentum component, p cos 6 Z , an observable on which the charge ratio is ex- 
pected to depend in a simple way [7]. 
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Table 1: Charge ratio R and relative statistical (stat.) and systematic (syst.) uncertainties in 
bins of p (in GeV/c), for surface data and both analyses of underground data. The relative 
uncertainties are expressed in %. 



V 


2006 surface 


2008 global-muon 


2008 standalone 


-muon 


range 


R 


stat. 


syst. 


R 


stat. syst. 


R 


stat. 


syst. 


5-10 


1.249 


2.3 


1.3 


— 


— — 


— 


— 


— 


10- 20 


1.279 


0.5 


1.5 


— 


— — 


— 


— 


— 


20- 30 


1.276 


0.7 


2.1 


— 


— — 


— 


— 


— 


30- 50 


1.279 


0.9 


2.6 


1.268 


1.2 2.1 


1.287 


0.5 


1.5 


50- 70 


1.285 


1.6 


3.4 


1.302 


1.2 0.6 


1.274 


0.5 


0.8 


70 - 100 


1.223 


2.1 


5.1 


1.274 


0.9 0.7 


1.272 


0.4 


0.9 


100 - 200 


1.287 


2.4 


8.9 


1.280 


0.8 0.3 


1.298 


0.3 


0.6 


200-400 


— 


— 


— 


1.295 


1.6 1.3 


1.305 


0.8 


1.4 


>400 


— 


— 


— 


1.349 


3.5 3.5 


1.350 


2.2 


6.0 



7 Results 

The results of the three analyses are shown in Fig. [61(a), as a function of the muon momentum. 
In the region where the results overlap, agreement between them is good, so the individual 
analyses are combined using a standard prescription [29 1. Within each analysis, some system- 
atic uncertainties are assumed to be correlated between momentum bins: trigger efficiency, 
momentum scale, charge misassignment and asymmetries in the muon losses due to the de- 
tector acceptance. In the global and standalone-muon analyses, systematic uncertainties from 
material densities, event selection, alignment, and magnetic field, are mostly uncorrelated be- 
tween momentum bins, and are treated as fully uncorrelated. On the other hand, they are 
correlated between the two analyses. 

The combined data points are given in Table [2] as a function of p and p cos 8 Z . They are shown 
in Fig. [61(a) as a function of p, and in Fig. |6](b) as a function of p cos 8 Z . 

Table 2: The muon charge ratio R from the combination of all three CMS analyses, as a func- 
tion of p and p cos 9 Z , in GeV/c, together with the combined statistical and systematic relative 
uncertainty, in %. 



p range 


(P) 


R 


uncertainty 


p cos 9 Z range 


(pcos0 z ) 


R 


uncertainty 


5-10 


7.0 


1.250 


2.45 


2.5- 10 


5.3 


1 27 A 


0.99 


10- 20 


13.7 


1.277 


0.85 


10- 20 


13.6 


1.251 


1.26 


20- 30 


24.2 


1.276 


1.34 


20- 30 


24.1 


1.262 


1.88 


30- 50 


37.8 


1.279 


1.10 


30- 50 


37.7 


1.292 


1.27 


50- 70 


58.5 


1.275 


0.54 


50- 70 


58.4 


1.267 


0.71 


70 - 100 


82.5 


1.275 


0.68 


70 - 100 


82.4 


1.289 


0.70 


100 - 200 


134.0 


1.292 


0.52 


100 - 200 


133.1 


1.292 


0.72 


200-400 


265.8 


1.308 


1.29 


200-400 


264.0 


1.330 


1.99 


>400 


698.0 


1.321 


3.98 


>400 


654.0 


1.378 


6.04 
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7 Results 
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(a) (b) 

Figure 6: (a) The three CMS results, and their combination, as a function of the muon mo- 
mentum. Data points are placed at the bin average, with the points from the standalone and 
global-muon analyses offset horizontally by ±10% for clarity, (b) The CMS result, as a function 
of the vertical component of the muon momentum, together with some previous measurements 
and a fit of the pion-kaon model to the CMS data. 

7.1 Charge ratio below 100 GeV/c 

In the region p < 100 GeV/c there are measurements in six p bins. Three bins are covered by 
all three analyses, with the surface-based MTCC analysis extending the reach to three lower- 
momentum bins. These twelve data points are combined into a single value of the charge ratio 
using the same prescription and scenario for correlations as for the overall combination de- 
scribed in the above section. This yields a charge ratio of 1.2766 ± 0.0032 (stat.) ± 0.0032 (syst), 
with a ^ 2 /ndf = 7.3/11, in good agreement with previous measurements |2]-|5l and represent- 
ing a significant improvement in precision. 

Repeating this fit in the pcos9 z region below 100 GeV/c yields a charge ratio of 1.2772 ± 
0.0032 (stat.) ± 0.0036 (syst.), with a ;t 2 /ndf = 15.3/11. The higher # 2 /ndf indicates that the 
data in this p cos 9 Z region have a lower probability of being consistent with a flat charge ra- 
tio. Fitting just the region p cos 9 Z < 70 GeV/c yields a charge ratio of 1.2728 ± 0.0039 (stat.) ± 
0.0040 (syst.) with a ^ 2 /ndf = 4.0/8, consistent with the flat charge ratio hypothesis. 

7.2 Charge ratio in the 5GeV/c to 1 TeV/c momentum range 

Considering the full p cos 9 Z range measured, a rise in the charge ratio is seen, as shown in 
Fig. [6] (b). Comparing to previous measurements in the same momentum ranges, the CMS 
results agree well where there is overlap: with the L3+C measurement [5] below 400 GeV/c, and 
with the UTAH [1], MINOS [6] and OPERA [8] measurements above 400 GeV/c. Measurements 
by other experiments in the range 5-20 GeV/c (2H3 133 are not shown in the plot; they are 
consistent with the constant value fitted in the CMS data. 



Models of cosmic ray showers provide an explanation for the rise in charge ratio at higher mo- 
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mentum. Based on the quark content of protons, and on the observation that primary cosmic 
ray particles are mostly positive, the ratio n + / n~ is predicted to be around 1.27 pl fl . Due to 
the phenomena of associated production, the charge ratio of strange particles such as kaons is 
expected to be even higher. 

The expected muon spectrum has been parametrized [32] based on the interactions of primary 
cosmic ray particles and on the decays of secondary particles, and from this parametrization, 
the charge ratio can be extracted as a function of the fractions of all pion and kaon decays 
that yield positive muons, f K and fa, respectively. These constants are not known a priori, and 
must be inferred from data. 

A fit performed to the combined CMS charge ratio measurement in the entire pcos9 z re- 
gion, with a fixed relative amount of kaon production [32|, yields f n = 0.553 ± 0.005, and 
fa = 0.66 ± 0.06, with a ^ 2 /ndf = 7.8/7. Figure |6J(b) shows the fit to CMS data only, together 
with a fit performed on some previous measurements by L3+C and MINOS 0. 

8 Conclusions 

We have measured the flux ratio of positive- to negative-charge cosmic ray muons, as a func- 
tion of the muon momentum and its vertical component, using data collected by the CMS 
experiment in 2006 and 2008. The result is in agreement with previous measurements by un- 
derground experiments. This is the most precise measurement of the charge ratio in the mo- 
mentum region below 0.5 TeV/c. It is also the first physics measurement using muons with the 
complete CMS detector. 
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